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We discuss functional-integral approaches to far-from-equilibrium quan- 
tum many-body dynamics. Specific techniques considered include the 
two-particle-irreducible efi^ective action and the real-time flow-equation 
approach. Different applications, including equilibration after a sudden 
parameter change and non-equilibrium critical phenomena, illustrate the 
potential of these methods. 

1. Introduction 

Time evolution of many-body systems in which effects of quantum physics 
play an important role belong to the least understood physical phenomena. 
Quantum effects prevent a fully deterministic description of time evolution 
and they render the problem mathematically difficult. Two main aspects 
constitute the complexity connected with non-equilibrium phenomena and 
their theoretical description: long-time evolution and strong correlations. 
These aspects can not be considered independently from each other. Due 
to the enormous progress in the field of ultracold atomic quantum gases 
in recent years, quantum many-body dynamics can now be very precisely 
controlled and probed and has become a timely topic for laboratory studies. 
On the theory side, the exponential growth of available computing power 
has brought quantum dynamics of many-body systems into the reach of 
large-scale calculations while sophisticated advanced analytical approaches 
open the way to a deeper understanding of quantum many-body evolution. 
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In situations close to equilibrium, one can treat non-equilibrium dy- 
namics as a perturbation of the equilibrium state. In different ways, most 
presently available field-theoretical methods rely to a certain extent on a 
perturbation expansion and either require a weak coupling or break down 
at large times of the evolution. This feature is inherent in principle, and it 
poses a serious problem when applying present results to quantum many- 
body dynamics far from equilibrium or to strongly coupled systems. The 
escape from this restriction is generally seen in refined approximations 
which allow to take complex correlations into account over a sufficiently 
long period of evolution before they cease haveing any further influence. 
Finding suitable approximations beyond leading-order perturbation theory 
can become technically involved. In this respect, functional-integral ap- 
proaches including the two-particle irreducible (2PI) effective-action and 
the real-time functional flow equation nicithods for quantum field dynam- 
ics represent successful ways of optimizing outcome versus effort. A cen- 
tral feature of these approaches is that the resulting dynamical equations 
conserve crucial quantities like energy, irrespective of the chosen approxi- 
mation. Most importantly, the methods allow approximations beyond an 
expansion in powers of the interaction strength. This review summarizes 
briefly the main aspects of the formalism and presents three examples for 
applications in which we compare with other methods and give speciflc 
predictions amenable to experimental investigation. 

2. Non-equilibrium quantum field theory 

Assuming basic knowledge about Feynman path integrals, we sketch the 
functional-integral formulation of real-time quantum fleld theory. To de- 
rive conserving many-body dynamic equations beyond mean-fleld order, we 
introduce the two-particle irreducible effective action as well as the real-time 
flow-equation approach. 

2.1. Schwinger-Keldysh closed time path 

We will consider initial-value problems, assuming that the many-body 

state is given by a density matrix /^(to) at some initial time t = Iq. 
The Schrodinger- and Heisenberg-picture evolutions of expectation val- 
ues are related to each other by the time evolution operator U{t,t') = 
Texp{— i J^, dt" H{t")} determined by the Hamiltonian H (we use h = 1), 

= Tr [p{t)Oit)] = Tr [p{to)U^it, to)0{t)U{t, to)] • (1) 



Functional-integral approaches to quantum many-body dynamics 



3 



CTP C: 



+ branch 



- branch 



Fig. 1. Schwinger-Keldysh closed time path C. The green dots indicate the times xq 
and yo for an example two-point function (^{x)^{y)) , see text. The branches are drawn 
above and below the time axis to make them separately visible. 



We use four- vector notation x = (xo,x) with time t = Xq and space coor- 
dinate X. The operator O is usually a product of field operators evaluated 
at different times such that, e.g., taking xq > yo, 

=Tr[p(to)C/^(xo,to)$(x){/(a;o,2/o)*(y)C/(yo,io)] ■ (2) 

The product of different time evolution operators and field operators eval- 
uates along a closed time path (CTP), Fig. [TJ If one chooses all operator 
times to lie on the -I- branch, the path integral gives the expectation value of 
the time-ordered product. In the path-integral approach, all evolution op- 
erators U are replaced by path integrals. Hence, the path integral consists 
of a product of integrals, one for each point along the CTP. 



2.2. Generating functional for correlation functions 

Within the Schwinger-Keldysh approach, the fundamental entity of non- 
equilibrium dynamics is the trace over the initial density operator multiplied 
by the forward and backward time evolution operators. Rewriting the latter 
as path integrals, the generating functional for correlation functions reads 

ZpJJ] =Tr \p{h)rce^fc^'^] ^ f 2?^p[^+, ^„ ]e«(5d^]+/c ^'.), (3) 



with p[ip'^,ipQ] — {'Po\p{to)\ipQ) ■ The external source field J{x) turns the 
path integral into a generating functional for correlation functions, similarly 
as in the (grand) canonical partition function in equilibrium physics, and 
7c denotes operator ordering along the CTP C. Furthermore, Jip = 
Jcd''+^xJ{x)ip{x) = /d'^x[//;""dxoJ+(2;)^+(x) -f //_^^^d:EoJ-(x)(p-(2;)], 
and Vip = nx['^V'^(2:)][d'/5"'"(a:)], with [dip^] = d(p^ {to , :x.) . \ipQ) are 
eigen states of ^{to, x), and Sc is a sum of two actions such that the overall 
time integral runs along the CTP. The generating functional allows, e.g., 
the field expectation value (/)(x) — {^{x)) to be written as 



' 6J{x) 



- Z^-i / ^(x) p[ip+ , ^0 ] . (4) 

j=o 
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Due to causality, the CTP extends only to the maximum time to be eval- 
uated in a particular n-point function. At later times, the sources can be 
set to zero such that the time evolution operators on the corresponding + 
and — branches cancel by unitarity. 



2.3. The two-particle irreducible (2PI) effective action 

The 2PI effective action is obtained by a double Legendre transform of 
—ihiZp^ with respect to J and to a further two-point source Rabix,y), 

r[<^,G] = r«[0]-i / Rab{x,y)[My)M^) + GUy,^)] (5) 

^ J xy 

where = dxo / d'^x, and r^[(^] = M^^[J] — / Ja4>a is the first Legen- 
dre transform of the Schwinger functional — — ilnZ^ for the mod- 
ified classical action Sq[(p] — Sc[(p] + j^ipRip/2. The classical source 
fields J and R act as Lagrange multipliers for the extremization of the 
2PI effective action with respect to 4> '"^^d the two-point Green function 
Gab{x,y) = {Tc'^a{x)^b{y)) - 4>a{x)(t)b{y), under the constraints Q and 
/ 5Ja{x)5Jb{y)\j=o = Gab{x,y). In turn, the action functional ^ al- 
lows to determine the correlation functions (j) and G from the Hamiltonian 
variational conditions 

— — — = -Ja(a;)- / Rab{x,y)(f>b{y), 77^—7 ^ ^ -7;Rab(x,y). (6) 

0(Pa[X) Jyfi dGab(x,y) 2 

Presuming a non- vanishing (/>, a saddle-point expansion can be used to 
write the 2PI effective action in terms of a one- loop part and a "rest" r2 , 

r[(/.,G] =5cM + ^Tr(lnG-i + GoiG) +r2[0,G]+const., (7) 

with free inverse propagator zGq ^^(x,?/) = S'^Sc[(j}]/S(l)a{x)S(l}biy)- For R = 
0, r2 yields the self-energy 

I],fc(x,y;(/.,G) = 2z^^[^ (8) 
dGab{x,y) 

that contains all information about scattering in the dynamics. The dy- 
namical evolution of G is governed by the Dyson equation obtained from 
the second equation in ([6|): 

G-\x,y)^Go\x,y)-J:{x,y). (9) 
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rj'|»»[(|),G] = 



r,""-" [G] : 




Fig. 2. (Left panel) Diagrammatic representation of the two- and three-loop diagrams 
contributing to r2 [(f>, G] in the 2PI effective action, for a system with elastic two-to- 
two collisions. The bare interaction vertices are drawn as black dots. (Right panel) 
Diagrammatic representation of the leading order (LO) and next-to-leading order (NLO) 
contributions in the l/A/'-expansion, to the 2PI part r2[4>, G] of the 2PI effective action, 
for = 0. The thick blue lines represent 2-point functions Gab{^,y) and the wiggly 
lines a single vertex channel. At each end of the wiggly lines, it is summed over the field 
indices of the G lines ending there, and at the whole vertex it is integrated over time 
and space. 



Upon multiplication by G this becomes a time evolution equation for the 
two-point Green function G, 

/ G(7L(a;, z)Gcb(2;,y) = 5a6'5c(a; - y) + / Y.ac{x, z)Gcb{z,y), (10) 

J z J z 

with 5c{x ~ y) — 5c{xq — yo)S{'x. — y). A corresponding equation can be 
obtained for ipa- Since the self-energy E is one-particle irreducible, and since 
taking the derivative with respect to G corresponds to opening a propagator 
line, it follows that T2 must consist of closed 2PI diagrams only, those which 
do not fall apart on opening two propagator lines, see Fig. This forms 
the central result that the 2PI effective action is given, besides the terms 
(O, by a series of all closed 2PI diagrams which can be formed from the full 
propagator G, the bare vertices defined by the classical action, and at most 
two external field insertions 0. We conclude that the 2PI effective action 
approach yields a closed set of dynamic equations for (p and G provided 
that the correlations in the initial state po can be fully encoded in these 
functions. This is the case for Gaussian initial conditions for which all n- 
point correlation functions with n > 3 can be expressed in terms of (f) and 
G0 Such higher-order correlations, however, build up during the evolution 
and are implicitly accounted for in the dynamic equations for and G. 

The 2PI effective action has been introduced to solid-state theory (there 
called ^-functional) [TJ [2] and later to relativistic quantum field theory 
[3]. See also Refs. [4j [5]. Applications to scalar relativistic as well as 

Tor non-Gaussian initial states a straightforward generalization of the approach involv- 
ing the nPI effective action is at hand |17| . 
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gauge theories can be found in Refs. [BHTT]. to non-relativistic systems, in 
particular to ultracold gases, in Refs. [T8U26) . 

2.4. Functional flow-equation approach 

We briefly sketch an alternative approach which is based on functional 
renormalization group (RG) techniques |27f[5^ . For more details on this 
approach, see Refs. [5S1|33]. Dynamic equations can be derived which are 
similar in structure to the equations obtained from the 2PI effective action. 
The key idea of the approach is to consider the generating functional for 
Green functions where all times are smaller than a maximum time t, im- 
plying the CTP to be cut off at r. The corresponding generating functional 
Zr is then defined in terms of the full generating functional Z — Zpg (j3]) as 



where the function Rt is chosen such that it suppresses the fields, i.e., 
d/SJa, for all times t > t. This does not fix Rt in a unique way, and one 
choice is —iRr^abi^, y) — for xq — i/q > t, x. = y, a ^ and zero 
otherwise, corresponding to Rt^oo = i? = in the 2PI approach. 

The time evolution of correlation functions is now derived from that of 
the Schwinger functional Wr = —ilnZT. It is more convenient, however, 
to consider the time evolution of the effective action 



An exact functional RG or fiow equation for the r-dependent effective action 
can be derived in the compact form 



where rV'^ = (5"r^/(50)". Eq. ^ is analo gous to the Wetterich functional 
fiow equation |27| used extensively with regulators in momentum and/or 
frequency space to describe strongly correlated systems near equilibrium 
P7H5^ . To obtain a practically solvable set of dynamic equations, one 
derives the flow equation for the proper n-point Green function Fr"'' by 
taking the nth field derivative of Eq. This scheme yields a set of 

coupled integro-differential equations for the Ft"\ and eventually for the 
connected n-point functions G^"-* , including the Schwinger-Dyson equation 
for the two-point Green function G. Its power lies in both a physically 




(11) 




(12) 




(13) 
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motivated evaluation of the effective action and the possibihty to derive 
non-perturbative equations of motion in a compact form [53l [54] . 

3. The 2PI effective action beyond mean-field order 

In order to evaluate the dynamic equation ()10|) , details about the self-energy 
E are required, and these are, in all interesting cases, only available to a 
certain approximation. The natural expansion of T2 is in terms of 2PI 
closed loop diagrams involving only bare vertices and full propagators G. 
This expansion can be truncated at any order, e.g., of powers of the bare 
coupling g or of the number of loops. 

3.1. Mean-field and quantum Boltzmann truncations 

Mean-field approximations that lead to the Gross-Pitaevskii and the 
Hartree-Fock-Bogoliubov (HFB) equations emerge as leading-order trun- 
cations of the 2PI effective action — 0, and T2 — double-bubble 
0{g) diagram, respectively, see Fig. [5]). Higher truncations take into ac- 
count non-mean-field effects of the collisional interactions. The closest con- 
nection to well-known formulations which account for scattering is obtained 
within the leading-order truncation in powers of the bare coupling g, taking 
into account the "Basketball" diagram in r2 . At this level of approximation 
the resulting dynamic equations can be reduced, after introducing quasi- 
particles at the mean-field level, and a Markov approximation within the 
scattering integral, to the well-known quantum Boltzmann equations for the 
single-particle densities, see, e.g., Ref. [21] ■ The Basketball diagram pic- 
torially illustrates that it describes two-to-two scattering with momentum 
conservation, and that it is of second order in the coupling. 

3.2. Conservation laws and long-time stability 

The most prominent advantage of the 2PI approach over, e.g., equations 
derived from the IPI effective action or the BBGKY hierarchy of coupled 
equations for cumulants, is that, whatever truncation is chosen, the re- 
sulting dynamic equations respect the conservation of energy and particle 
number. This property is due to the self-consistent determination of G, 
and (j) in the case of bosons. As a result, provided an implementation with 
sufficient numerical precision, the equations do not lead to a secular evolu- 
tion or the emergence of negative occupation numbers. This is particularly 
important for approximations which include the effect of scattering and 
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allow to describe long-time evolution including equilibration and thermal- 
ization [14) . 

3.3. Non-perturbative expansions beyond mean-field order 

While 0{g") truncations of r2 render the approach useful for perturbative 
approximations beyond mean- field-order, resummations of infinite classes 
of diagrams are possible which extend the applicability beyond the cou- 
pling expansion. In the following, we will review a few applications of a 
resummation procedure corresponding to the next-to-leading order (NLO) 
truncation of an expansion of r2 in inverse powers of the number of in- 
ternal field degrees of freedom TV [H [7], see Fig. [2] and Refs. (HHIIl EOHSS] 
for non-equilibrium applications. This procedure takes into account clas- 
sical statistical fluctuations to infinite order in the coupling and includes 
quantum fluctuations to leading order in the quantum contribution to the 
interaction vertex |23j . Hence, in this truncation, the dynamic equations 
are particularly useful when studying the dynamics of systems in which 
quantum statistical fluctuations are suppressed. This is, e.g., generically 
the case in Bose gases at low energies in which the bulk of particles occu- 
pies a few excited modes near zero energy. In the context of Kadanoff-Baym 
equations the NLO l/Af expansion is known as GW-approximation [551 [5S] 
and has been the subject of non-equilibrium studies recently [57H59j . 

4. Applications 

We present three examples where the dynamic equations derived from a 
non-perturbatively approximated 2PI effective action provide insight into 
the many-body time evolution beyond the mean-field as well as leading- 
order quantum-Boltzmann equations. 

4.1. ID lattice gas: 2PI versus exact dynamics 

In the first example, we consider a coupled few-mode Bose system as, 
e.g., a one-dimensional (ID) lattice gas in the tight-binding approxima- 
tion. Its full quantum dynamics can be computed exactly or simulated 
in the quasiclassical regime, and we compare such results with predictions 
obtained within the 2PI approach. The system is defined by the Bose- 
Hubbard Lagrangian C(n,t) ^ [|^';(t)9t*„(t) + J^-* (t)5'„+i(t) + c.c] - 
e„5'*(t)5'„(t) - ^^['^*^{t)^ n{t)f , n 1, A^s being the lattice site index. 
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Fig. 3. (a) Time evolution of the Josephson contact between Ns = 2 bosonic modes. 
The three columns show the same quantities for a different total number of atoms 
N = 20, 40, and 80, respectively. NU/J = 4 throughout. Initially, all atoms are in 
a Bose-Einstein condensate in mode 1. First line: Total number of atoms A'^tot.i in 
mode 1. Due to number conservation, A'^tot,2 = N — A^tot,!- Second and third lines: 
Condensate fractions Nc,„/N at sites n = 1,2. Fourth line: Nc/N = JZ^Nc,n/N . 
expansion to NLO: Thick solid line. Blue dashed lines: Exact calculation. Green dash- 
dotted curves: l/J\f expansion in 2nd-order coupling approximation. Red dotted curves: 
Hartree-Fock-Bogoliubov approximation, (b) Upper: Time evolution similar to (a) for 
different numbers of sites Ns = 2 (black), 3 (red), 5 (blue) and 10 (green). The other 
parameters are chosen as NU/NsJ = 5 and N/Ns = 10. Solid line: l/J\f expansion 
to NLO. Dashed line: quasi-exact semi-classical simulations. Lower: Evolution of the 
condensate fraction and of the total population (inset) per site for A^^ = 2. 



For Ns = 2, the system of N particles is equivalent to a Josephson junc- 
tion with coupling energy i?j = 2JN and charging energy Ec — 2U, cf., 
e.g., Refs. [60l [61]. In order to probe the accuracy of the NLO 1/J\f ap- 
proximation, we chose NU/ J to be larger than 1 and compare our results 
with the results of an exact numerical calculation. In Fig. [3^, the 2PI 
evolution in NLO 1 /TV approximation is compared to the HFB and second- 
order coupling approximations, see Ref. |22j for further details. The time 
evolution shows different characteristic periods. At early times, the con- 
densate Rabi oscillates coherently between the lattice wells with frequency 
= 2E}/N = 4 J. Only a small number of atoms is scattered from the con- 
densate fraction into excited modes. Lateron, atoms are exchanged between 
the condensate and the non-condensate modes of the gas. These processes 
lead to a rapid destruction of the condensate fraction and to damping of 
the Rabi oscillations. 

In Fig.[5}D, we compare the evolution of the total condensate fraction ac- 
cording to the 2PI NLO 1/Af equations with the quasi-exact semi-classical 
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simulations ("Truncated Wigner Approximation"), for Ng up to 10. These 
results illustrate that the field-theoretic non-perturbative resummation is 
capable of qualitatively describing the evolution. Note that, in contrast to 
this, coupling expansions generically break down at large times. Quantita- 
tive agreement converges at most slowly with increasing number of sites. 
We furthermore find that the total number of particles at a particular site is 
described much better, by the NLO equations, than the condensate fraction 
alone, which is strongly overestimated, see Fig. for an example. This 
indicates that, in the NLO l/Af approximation, the phase spreading in the 
different modes is captured only to a limited extent. 

4.2. Non-thermal equilibration of a ID Fermi gas 

As a second example, we study the long-time evolution and equili- 
bration of a ID Fermi gas [621 |63] . containing two spin components 
a G {tii} that mutually interact through local repulsive s-wave colli- 
sions described by the Hamiltonian H = f dx['i/l^{x){—d^/2m)'ia{x) + 
igai3/2)^i{x)^l{x)'i'p{x)'^a{x)l with g„/3 = (1 - S^p) 'inaiB/m, am be- 
ing the ID scattering length. This model is considered integrable in the 
sense that it has as many conserved quantities as there are degrees of free- 
dom [64]. As a consequence, the gas is expected not to be described by a 
grand-canonical density matrix at large times |65fl67j . 

An extensive discussion has recently focused on the question under 
which circumstances a generalized Gibbs ensemble (GGE) can describe the 
long-time properties of, in particular, integrable systems p5 l [68H72] . From 
the experimental point of view, it appears difficult to decide on which gen- 
eral grounds one may be in the position to falsify the equilibration to a 
particular final density matrix. We will show that, instead of restricting 
a measurement to a single point in time, it is useful to look for a possible 
deviation from the fiuctuation-dissipation theorem. The results indicate 
that a degenerate ID Fermi gas defined by the Hamiltonian H approaches 
a non-thermal state at large times. 

In Fig. HJj, the time-evolution of six single-particle momentum modes 
iT'titj |Pi|) = n-iity |Pi|) shown, for an unpolarized gas following an in- 
teraction quench at to from 7 ~ mg/niD = to 7 = 4, for different initial 
momentum distributions n^(|p|) — nf{to,\p\) (Runs A and B in Fig. 
The equilibration process shown in Fig. |3)d is characterized by a short-time 
scale of dephasing depending on the width of the initial momentum distri- 
bution and a long-time scale determined by scattering. The different initial 
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Fig. 4. Equilibration of a one-dimensional Fermi gas starting from different initial 
momentum distributions, (a) Initial (t = 0) and final (t = 10 mnj"^) distributions 
^fi'tj IpI) = "■4,(i) IpI) in runs A and B have the same total particle number and energy. 
In runs C and D, the particle number is the same while the energies are lower than 
in B. (b) Occupation numbers n^(4, |p|) as a function of time t for momentum modes 
Pi = s'm[iiT/Ns]/as, Ns = 128. (c) Inverse-slope functions a = ln[l/n.|-(p) — 1] for runs 
B, C, and D at late times. The black dashed lines are Fermi-Dirac fits to the lowest 14 
momentum modes from which the temperatures and chemical potentials are extracted. 
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Fig. 5. 
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(a) Late-time momentum distributions equivalent to those shown in Fig. [J};, (b) 



in runs C and D for five of the momentum modes pi, see legend in lower right panel. 
Lower row: Inverse-slope function of fractions / of the statistical correlation function 
F divided by the spectral function p at Xq = IS.Qn^j^m for the same five momentum 
modes, see main text. Black lines indicate Fermi-Dirac distributions with the same /3 
and /i as in Fig. |4j;. In run D, the system does not thermalize. 



distributions of runs A and B have the same particle number and energy 
such that both runs approach the same final state. The initial momentum 
distributions in runs C and D further shown in Fig.|3^ contain the same par- 
ticle number as runs A and B but have lower energy. We find that also runs 
C and D reach stationary momentum distributions which, however, are no 
longer given by a Fermi-Dirac function over the entire range of momenta. 
As seen in Fig. |4]:, the lower momentum modes appear thermalized while 
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the higher momenta in runs C and D remain overpopulated as compared 
to the exponential fall-off of the Fermi-Dirac distribution. 

Fig. [5^ shows the n-[.(|p|) for runs C and D on a double- logarithmic scale. 
The high-momentum tails are characterized by a power-law n(p) oc p^** 
with K ~ 4.4. This exponent does not change when we consider trun- 
cations of the full 2PI loop expansion that contain all diagrams up to 
order or g^. This non- Fermi-Dirac momentum distribution in itself 
is not a proof for non-thermal equilibration. However, the final state is 
found to also violate the fluctuation dissipation theorem which points to 
a non- (grand-) canonical character. Given a grand-canonical many-body 
state p — Z^^ exp[(3{H — pN)] where N is the total-particle-number op- 
erator, the fluctuation-dissipation theorem states that the statistical cor- 
relation function FaaiXo;uj,p) — J ds e:sq){iu!s)Faa{XQ + s/2,Xo — s/2; 
p), where Faa{t,t';p) = ([4'J,(t, p), ^'^(t', p)])/2, and the spectral function 
Paait,t';p) = (t,p),^'„(i',p)}) are related by 

Fq.q(Xo; w,p) = "«[l/2 - riFD(i^ - ^J.)]paa{Xo;^^,v)■ (14) 

In Fig. O we depict, at the late time Xq = 18.9 n^^m, the fraction 
/ = i_F|-|-(Xo; Lu, p)/p^^(Xo; w, p) as a function of the frequency uj for five dif- 
ferent momentum modes p^. The lower left (right) panel shows the inverse- 
slope function ln[(l/2 — f)^^ — 1] of / for run C (D). / is shown in a region 
including the spectral peaks where the argument of the logarithm is posi- 
tive. In run C, over the region of relevant to, this function is a straight line 
corresponding to a Fermi-Dirac function /. Hence, the system is thermal- 
ized over the depicted range of energies, despite the power-law tail in run 
C found for Wp > 30 wf, see Fig. St. This can be understood by consider- 
ing the spectral function in Fig. [5]d (upper row) on a logarithmic scale. A 
second peak at negative frequencies, not present in an ideal gas, picks up 
extra contributions from the Fermi sea. It is this second peak which causes 
the power-law overpopulation at high momenta, in the same way as of the 
BCS zero-temperature depletion of the Fermi sea in a weakly interacting 
system. Hence, in run C, the system thermalizes to a state within a grand- 
canonical ensemble, with the eigenmodes of the strongly interacting system 
being superpositions of particles and holes. Run D, however, performed 
at even lower energy, shows that the system does in general not thermal- 
ize to a grand-canonical ensemble. As shown in the lower right panel of 
Fig. the result of this run violates the fluctuation-dissipation theorem. 
Although the momentum overpopulation is again largely produced by the 
contributions from the Fermi sea, also the fraction / shows a power-law 
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tail ^ p~^- In conclusion, the 2PI equations describe non-thermal equili- 
bration of a ID Fermi gas with positive contact interactions, at sufficiently 
low total energies, to a state violating the fluctuation-dissipation theorem 
for a grand-canonical ensemble. 

4.3. Non-thermal fixed points and quantum turbulence 

In the last example, we discuss non-equilibrium fixed points of a Bose gas. 
The dynamic equations derived from the non-perturbatively approximated 
2PI effective action can be analyzed with respect to stationary scaling solu- 
tions. The idea is to ask for a solution of the dynamic equation (jlOp that is 
translation-invariant in time and space, and is a scaling function in the fol- 
lowing way: The frequency-momentum dependent statistical function (see 
previous section) obeys F{s''qQ, sq) — s~^~''F((7o, q), while the spectral 
function satisfies p{s^qo, sq) = s~^p{qo,q), the latter in accordance with 
the commutation relations. The scaling law for F implies the single-particle 
spectrum n(k) to scale as 

n{sk) = s~^n{k), (15) 

where s is some positive real number and ( = k -\- 2 ~ z. Within kinetic 
theory, where the quantum Boltzmann equation takes the role of the equa- 
tion of motion for n(k, t), scaling laws can be derived by use of Zacharov 
integral transformations, 'little miracles of wave-turbulence theory' |73) . 
These transformations allow to rewrite the complicated scattering integral 
in such a way that scaling exponents of solutions other than the Rayleigh- 
Jeans distribution n ~ T/fc^ can be read off. 

To determine the positive exponent ^ in the infrared (IR) regime, where 
n(k) ^ |k|^^ is large, an approach beyond kinetic theory is required. Recall 
that the 'Basketball' diagram in the 2nd line of Fig. [2] (left panel) gives the 
quantum Boltzmann equation. This diagram is proportional to a lower 
power of the occupation number n than all higher-order ones. Hence, the 
approximation becomes unreliable for n 1. In contrast, the NLO l/Af 
approximation or s-channel resummation involves diagrams to all orders in 
the coupling and in n. While, in the regime of large wave numbers, the 
approach goes over into the kinetic description of 'weak wave turbulence', 
an effectively renormalized many-body T-matrix modifies the scaling law 
in the IR. Details of this and of the derivation of the scaling exponent C are 
given in Ref. 25 . 

The renormalization of the coupling is, in the NLO l/7V-approximation, 
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Fig. 6. Left panel: Phase angle (color scale) at four times during a single run of the 
simulations in 2 dimensions. The white spots mark vortex cores where the density falls 
below 5% of the mean density. Shown times are (in lattice units): 1. I = 26 (top left): 
Ordered phase shortly after initial preparation. 2. t = 820 (top right): After creation of 
vortex-antivortex pairs. 3. t = 6550 (bottom left): During critical slowing down of the 
vortex-antivortex annihilation. 4. t = 10^ (bottom right): Low-density vortex-antivortex 
pairs before final thermalization. Right panel: Mode occupation numbers as functions 
of the radial momentum k, for the four diflperent times of the run in 2 dimensions shown 
on the left. Note the double-logarithmic scale. An initial broadening of the momentum 
distribution is followed by the early development of a scaling n(fe) ~ k~^ and the later 
emergence of the bimodal scaling with n{k) ~ k~^ at larger wave numbers. Scaling laws 
are indicated by thin lines. 

a consequence of the resummation of an infinite number of Feynman dia- 
grams contributing to the 2PI effective action in terms of a geometric se- 
ries [m [25l [74] . Physically, the renormalized T- matrix implies a reduction 
of the effective interaction strength in the IR regime of strongly occupied 
bosonic field modes [16]. As a consequence of the reduced interactions, 
single-particle occupation numbers rise towards smaller wave numbers in 
a steeper way than in the (ultraviolet) kinetic regime. The following IR 
scaling exponents were predicted in Ref. [25] : 

C = Cq =rf + 2, (16) 
( = (IR = d + 2 + z (17) 

where d is the number of spatial dimensions, cj ^ fc^, and Q (P) indicates 
that the quasiparticle number (energy) flows at constant rate in the cascade. 

Recently, these scaling exponents where shown, by means of simulations 
of the Gross-Pitaevskii equation, to be closely related to the phenomenon 
of quantum turbulence [75] . Given a generic out-of-equilibrium initial con- 
dition for an ultracold quantum gas, topological features such as quantized 
vortices or vortex lines may appear after a short evolution time. They are 
usually produced in great numbers, in 2D together with antivortices, to 
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comply with the overah angular momentum of the system, in 3D in the 
form of closed lines. While diminishing again in number due to pair annihi- 
lation or reconnection, clear power laws appear in the occupation number 
distribution, confirming the predictions obtained from the non-perturbative 
2PI approach, see Fig. |6l See Ref. [76] for a large- A/" relativistic study. 

5. Summary 

We have briefly reviewed functional-integral approaches to non-equilibrium 
quantum many-body dynamics, including the 2PI-effective action and the 
functional flow-equation method. These techniques provide a powerful and 
economic access to dynamics of strongly correlated systems, both analyt- 
ically and numerically. Approximations far beyond the mean-field level 
conserve crucial quantities like energy and particle number and allow to 
describe long-time dynamics and equilibration. With a few example ap- 
plications we illustrated the potential of the approach, demonstrating the 
description of equilibration both to thermal and non-thermal states, the 
evolution of uniform as well as inhomogeneous systems, of Bosons and 
Fermions. Also specific properties of non-equilibrium stationary states like 
critical scaling can be described, and the comparison with alternative meth- 
ods proves that the methods are useful for studying non-trivial topological 
solutions. 
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